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Abstract 

We study non-linear effects of radiative viscosity of npe matter in neutron 
stars for both direct Urea process and modified Urea process, and find that 
non-linear effects will decrease the ratio of radiative viscosity to bulk viscosity 
from 1.5 to 0.5 (for direct Urea process) and 0.375 (for modified Urea process). 
Which means that for small oscillations of neutron star, the large fraction 
of oscillation energy is emitted as neutrinos; but for large enough ones, bulk 
viscous dissipation dominates. 

Key words: neutron stars, bulk viscosity, radiative viscosity 
PACS: 97.60.Jd, 21.65.-f, 95.30.Cq 

1. Introduction 

Bulk viscosity can damp the density oscillations in compact stars, which 
could be excited at the time of their formation from supernovae explosions, 
or during the phase transitions Jl(, or due to instabilities result from the 
emission of gravitational waves 0, 3, 0, 0, 0] . As one of the most impor- 
tant transport coefficients, bulk viscosities of simple npe matter, of hyperon 
matter and even of quark matter, both in normal and superfuild states, have 
been extensively studied^, 3, 3, 03, EH, 03, 03, 03, 03, 03, 03, 03, 03, 0, S 



for more references see [221 ] . However, until recently, Sa'd et al. [23] demon- 
strated that there exists a new mechanism for dissipating the energy of stellar 
oscillations. They indicated that the mechanical energy of density pertur- 
bations is not only disspiated to heat via bulk viscosity, but also is radiated 
away via neutrinos. They named this new mechanism the radiative viscosity, 
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and found it is 1.5 times larger than the bulk viscosity to all Urea processes, 
both in nuclear matter and quark matter. 

The newly realized radiative viscosity was calculated only in the lowest 
order of b[ijT jijjj], which corresponds to the linear approximation. However, 
the density oscillations during the stellar evolution may arise to sufficiently 
large amplitude so that no-linear effects can no longer be neglected. We aim 
to study the non-linear effects of radiative viscosity in this paper. 

We will study non-linear effects of radiative viscosity of simple npe matter, 
both for direct Urea process (n — > p + e + V e , p + e — > n + u e ) and for modified 
Urea process (n+N — > p+N+e+V e , p+N+e — > n+N+u e ). As indicated by 
Lattimer et al. [24j, if the proton and electron Fermi momenta are too small 
compared with neutron Fermi momenta, the nucleon direct Urea process 
is forbidden because it is impossible to satisfy conservation of momentum. 
Under typical conditions, one finds that the ratio of the number density of 
protons to that of nucleons must exceed 1/9 for the process to be allowed. 

We not only make a numerical solution of radiative viscosity bue also 
study the ratio of the viscosity coefficient to the bulk one. By doing so, we 
hope to know which mechanism is more important under special condition. 

This paper is arranged as follows. In Sect. 2, we derive the formulae of 
radiative viscosity including non-linear effects. And radiative viscosity for a 
specific model of neutron star (NS) matter is calculated in Sect. 3. Finally, 
Sect. 4 presents the conclusions and discussions. 



2. The formulae of radiative viscosity 



In this section we derive the expression for radiative viscosity of neutron 
stars whose cores consist of simple npe matter using an approach similar 
to that of bulk viscosity calculations has been done by Wang and Lu 
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Sawyer [12L Madsen [13j and Gupta et al. [14 



We first recall how to solve chemical potential perturbations with stellar 
oscillations. Assume that the volume per unit mass, v, changes periodically 
in time according to the relation 



v(t) = vq + A-usin ( j = v + Sv(t), 



(1) 



where vq is the equilibrium volume, At; is the perturbation amplitude, and 
t is the period. 
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Taking the chemical potential difference as 



(2) 



where p n , p p and p e are the chemical potentials of the neutrons, protons and 
electrons. 5p can be expanded near the equilibrium according to 
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(note that 5p — in equilibrium) rij are particle numbers per unit mass and 



5n p = —5n n = 5n e 



(dn p /dt)dt. 



Using the thermodynamical relations 
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(5) 



and considering p { = rii/v (pi is the particle numbers per unit volume), one 
gets 



5a(t) = -A— sin f—\ + B [ (dp„/dt)dt, 

VO \ T J J 
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and P(t) is the pressure, and the net reaction rate is 



(6) 

(7) 
(8) 
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According to 25j, |26( and introducing z = Sp/nT, for direct Urea process 

7Uz 2 + 42(k 4 + 42z 6 



T d (T,5p)5p = e d (T,0)- 



457 



(10) 
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/ 1/3 

e d (T,0) = 3.3xl(T 14 p^ T 6 MeV 5 , (11) 



Po 



and for modified Urea process 



14, 680z 2 + 756(k 4 + 84(k 6 + 24z 8 

r m (T,<J/x)5/x = e m (T,0) , (12) 

11, old 

/ \ 1/3 

e m (T, 0) = 3.6 x 1(T 18 f T 8 MeV 5 , (13) 

where x p = p p / p^, is the ratio of the number density of protons to that of 
nucleons, and po = 0.16/m^ 3 . 

Given a specific equation of state (EOS) of npe matter, 5{i(t) can be 
calculated from eqs. (jSJ), ©, (|EJ) and (jUJ) numerically, which will be done 
in section 3. Once the time dependence of chemical potential difference is 
known, one can easily calculate the bulk and radiative viscous coefficient. 

For a periodic process, the expansion and contraction of the system will 
induce not only the dissipation of oscillation energy to heat, but also the loss 
of oscillation energy through an increasing of the neutrino emissivity. Bulk 
viscous coefficient ( and radiative viscous coefficient 1Z can be defined for the 
description of these dissipation mechanisms, respectively 23 



(4ss) = -- / dt(V-v) z , (14) 



T 



(4»> = - / dt(V-v)\ (15) 



r 



where v is the hydrodynamic velocity associated with the density oscillations. 
Using the continuity equation, one obtains 

C=-2(^)(|l) 2 (^) 2 , (16) 



K = 2 <M^) (t) • (17) 

where the energy dissipation is 

(4 SS ) = - / T(T,5fj.)5ndt, (18) 
Jo 
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and the neutrino emissivity caused by the oscillation of 8fi is 

(£ioss>= / [e(T,5n)-e(T,0)]dt, (19) 
Jo 



for direct Urea process [25 



e d (T,8^) =e d (T,0) I 1+ — 1, (20) 

and for modified Urea process 

e m (T, 5//) = e m (T, 0) I 1 + — — j. (21) 

3. Radiative viscosity for a specific model of NS matter 

3.1. EOS model 

For illustration, we use a phenomenological EOS proposed by Prakash et 



al. [27|. According to these authors, the nuclear energy is presented in the 
form 

E N ( Pb ,x p ) = E m ( Pb ) + S( Pb )(l - 2x p ) 2 , (22) 

where E^(p b , x p ) = Eiy (p b )(p b , x p = 1/2) is the energy of symmetric nuclear 
matter and S(p b ) is the symmetry energy Supposing the electron energy is 
E e (p b ,x p ), one has d[E^{p b ,x p ) + E e (p b ,x p )]/dx p — in (3 equilibrium. 

At the saturation density po = 0.16/m -3 the symmetry energy is mea- 
sured in laboratory, S(po) = 30MeV; while at higher p b it is still unknown. 



Prakash et al. 27] displayed S(p b ) in the following form 



S(p b ) = 13MeV [u 2/3 - F(u)] + S(p )F(u), (23) 

where u = p b /po and F{u) satifies the condition F(l) = 1. We employ 
their model I (F(u) = u) with the compression modulus of saturated nuclear 
matter K = 240MeV. This is a moderately stiff EOS, the maximum stellar 
mass for it is M — 1.977M Q . The direct Urea process is allowed only if 
x p > 1/9; in our EOS model it demands p b > 0.434fm -3 , which means the 
direct Urea process is completely forbidden in NS when its mass is smaller 
than M D = 1.358M . 
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Considering p p = p e = p b x p and p n = pb(l — x p ), one has 

(24) 





( dp' 




K dx P, 




f d5p 


pb 


\dx p 



B =-\if) ■ ( 25 ) 

b P / 

Thus, 5p(t) can be solved numerically, and the bulk and radiative viscosity 
can be calculated from eqs. ( fTBl) and ( fTTl) . respectively. 

3.2. The results 

The solid line in Fig.l and Fig. 2 shows radiative viscosity as function of 
relative volume perturbation amplitude for p& = 0.6/m -3 and p& = 0.3/m" 3 , 
note that the onset of the direct Urea process is p b = 0.434. As indicated 



in ref. [14j, the flat branches in Fig.l for the direct Urea process follows 
the T 4 law: when T increases by four orders (form 10~ 4 MeV to lMeV), vis- 
cosity increases by about 16 orders. In contrast, the flat branches in Fig. 2 
for the modified Urea process follows the T 6 law, and the radiative viscosity 
of modified Urea process is much lower than direct Urea process. However, 
the behavior of the onset of non-linear effect is similar for the two different 
processes: the lower the temperature, the smaller the relative volume pertur- 
bation amplitude is for non-linear effects to dominate. In the other words, 
as T decreases, non-linear effect becomes more and more important. Further 
more, we also present the results up to z 4 . As we all know, when calculated 
up to z 2 , one can only get the linear results; and non-linear effects are usually 
introduced by calculating up to z 4 . But by comparing solid and dash curves 
in Fig.l and Fig. 2, we can see that it's not enough just including the lowest 
order of non-linear effects. 

Fig. 3 presents the ratio of radiative viscosity to bulk viscosity as function 
of relative volume perturbation amplitude. For both direct Urea and modi- 
fied Urea, 71/ ( = 1.5 when Av/vq is sufficiently small, which is in agreement 
with the results given by Sa'd and Schaffher-Bielich [23j. And for both pro- 
cesses, non-linear effects lead to the decreasing of 7Z/(, and the lower the 
temperature, the more remarkable the non-linear effect is. Nevertheless, for 
direct Urea process, TZ/C can reach 0.5 when Av/v is large enough; but 
TZ/C can be equivalent to 0.375 for modified Urea process. Comparing the 
formulae (TTU|) with (1201) . and f)12p with (I2ip . one immediately know that the 
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ratios only depend on the average value of the net reaction rate and the in- 
crements of neutrino emissivity due to off-equilibrium, and EOS can't effect 
the results at all. It is also interesting to plot TZ/( as function of A/x/T 
(Fig. 4), where A/x is the perturbation amplitude of chemical potential dif- 
ference 5fi(t). It's undoubtable that the curves are overlapped for different 
values of temperature in Fig.4. 

4. Conclusions and discussions 

We have studied the non-linear effects of radiative viscosity of simple npe 
matter, both for direct Urea process and for modified Urea process, and found 
that non-linear effects decrease TZ/( from 1.5 (the linear scenario) to 0.5 (for 
direct Urea process) and 0.375 (for modified Urea process). Which means 
that in the linear scenario, only | times of the total oscillation energy is 
converted to heat and the large faction of the energy is emitted as neutrinos; 
in contrast, if the amplitude of the oscillation is large enough and the non- 
linear effects cannot be ignored, the main part of the total oscillation energy 
is converted to heat ( | for direct Urea process and ^- for modified process). 
Although the numerical results of the coefficients are calculated for a specific 
EOS, the ratio of the two types of viscosity coefficient is EOS-independent. 

In the case of superfluid npe matter, since superfluidity causes different 
effects to the net reaction rate and the increments of neutrino emissivity 
due to off-equilibrium, we expect that both in linear regime and in non- 
linear regime, the value of TZ/( will have great differences comparing with 
the normal EOS. This is the further work we set about. 
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Figure 1: Radiative viscosity as function of relative volume perturbation amplitude for 
t = 10 _3 s and pb = 0.6/to~ 3 , where the direct Urea process occurs. The dash curves are 
the results up to z 4 and the solid curves up to z 6 . The temperatures are 1CU 4 , 10~ 3 , 10~ 2 , 
10 _1 , and 1 MeV from bottom to top, respectively. 
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Figure 2: Radiative viscosity as function of relative volume perturbation amplitude for 
t = 10~ 3 s and pi, = 0.3/m~ 3 , where the direct Urea process is forbidden. The dash curves 
are the results up to z A and the solid curves up to z 8 . The temperatures are 10 -4 , 10 -3 , 
1(T 2 , 10 _1 , and 1 MeV from bottom to top. 
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Figure 3: lZ/( as function of relative volume perturbation amplitude for r = 10 3 s , 
Pb = 0.6/m~ 3 (solid curves) and pb = 0.3/m -3 (dash curves). The temperatures are 
1CT 4 , 1(T 3 , 1CT 2 , 10" 1 , and 1 McV from left to right. 
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Figure 4: lZ/( as function of lg(A/z/T) for r = 10~ 3 s , pb — 0.6/m -3 (solid curves) 
and pb = 0.3.fm~ 3 (dash curves), where Ap is the perturbation amplitude of chemical 
potential difference Sp(t). 
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